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Abstract 



In this paper some properties of decomposable semigroups and their 
associated ideals are studied. We give an efficient method to determine the 
decomposition of this kind of semigroups that improves the computation 
of some objects related with them. The concept of decomposable variety 
is introduced and a combinatorial characterization of decomposable semi- 
groups is proven. These improvements are applied to some semigroups 
studied in Algebraic Statistics. 

Introduction 

Let S C Z™ X Z/ciZ X ••• X Ijch^ be the semigroup generated by ^ = 
{ai, . . . , a„}, with positive integers ci, . . . , c/i, and such that SC\{—S) = {0}. Let 
k be a field and k[Xi, . . . , X„] the polynomial ring in n indeterminates. There 
exists a 5— grading of the ring k[Xi, . . . ,X„] assigning the S'-degree to the 
indeterminate Xi. Thus, the S'-degree of X" = X°^ ■ ■ ■ X^" is X]r=i "^i^* ^ 
we denote it by S'-degree(X"). The finitely generated S'-homogeneous binomial 
ideal generated by the set 



is called the ideal of semigroup S, denoted by Is (see [E]). In this paper we 
only consider binomial systems of generators. 

The ideal Is is equivalent to the congruence determined by 
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where p* : N" — > S* is the epimorphism p*{ai, . . . , a„) = Y^^=i '^i'^i ^ this 
case S is isomorphic to N"/p (see [131 Page 176]). In [5], [53] and [55] appear 
different techniques to compute (minimal) system of generators of semigroup 
ideals. 

A semigroup S is decomposable if it is the direct sum of some proper sub- 
semigroups, Si, such that 

Is=Is^+Is2+--- + Is,- (1) 

If a semigroup is not decomposable we say that the semigroup is irreducible. In 
this paper we assume that each semigroup Si is irreducible. 

If there exists a system of generators of Is that can be expressed as Bi U B2 
with Bi n B2 = 9 and the binomials in Bi have not indeterminates in common 
with the binomials in B2, then one has a decomposition of S (not necessarily 
irreducible) as described in ((U . One of the contributions of this paper is to give 
a method to obtain the decomposition into irreducibles of a semigroup without 
computing the ideal Is or a system of generators verifying the properties of Bi 
and i?2- 

Since the semigroup S has the property S Ci {~S) = {0}, by Nakayama's 
lemma (see [5]), all of the minimal generating sets of Is have the same cardi- 
nality. Denote by Betti{S) the set of 5— degrees of the elements of a minimal 
system of generators of Is- It is known that if a semigroup S is cancellative, 
it verifies the property S fl {—S) — {0} and it verifies the condition of ascend- 
ing chain (see [3]), then all the minimal systems of generators of Is have the 
same Betti{S) and the same number of binomials of 5'-degree m & S for all 
m S Betti{S). Every minimal system of generators of the ideal is called a 
Markov basis. We assume the semigroups in this work verify these conditions. 

The complexity for computing a semigroup ideal is simply exponential in 
the number of variables (see [22]) ■ The following rule is based in this setting: 
it is better to compute several ideals in fewer indeterminates that one in more 
indeterminates (see [26] and the references therein). In our case to compute Is 
we calculate the ideals /g. which are ideals in fewer indeterminates than Is- 

One of the most important objects in Algebraic Statistics (see [5] , [7] , [12] and 
[l9l about Algebraic Statistics) is the ideal associated to semigroups related with 
statistical models. In this context, one can find many examples of decomposable 
semigroups. Note that in the database ,141, 28% of the reducible graphical 
models are associated a decomposable semigroup. This justifies our study of 
decomposable semigroups. The uniqueness of the Markov bases (see [T] and 
[24) 1 is also a problem studied in Algebraic Statistics. 

Another application of this splitting in irreducible semigroups appears in 
the scope of non-unique factorizations invariants on monoids. The catenary 
and maximum of the Delta sets (see [TT] for the definition and properties of 
these invariants) is reached in the Betti elements of the semigroup (see [3] and 
the references therein). Since the set of Betti elements of the semigroup is the 
union of sets of the Betti elements of its irreducible factors, the study of these 
invariants reduces to the study of these invariants in simpler monoids. The same 
holds for the tame degree and the elasticity (see [TT] for its definition and prop- 
erties), because clearly from the definition of decomposition into irreducibles, 
the splitting also applies to the primitive elements of p (see [TBI'). 

Now, we summarize the content of this paper. 
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In section[TJ we introduce definitions, notations and some known results that 
are needed for understanding this work. 

In section [21 we show how the Hermite normal form can be used to obtain a 
diagonalization of a matrix that we call HNF-diagonal matrix. Algorithm [7] is 
the key that allows to compute easily the decomposition of a semigroup. 

In section[3J we characterize decomposable semigroups by means of the HNF- 
diagonalization of a matrix associated to S. This characterization allows us to 
study some properties of S by using the semigroups of its decomposition. Fur- 
thermore, we can also obtain semigroups that are isomorphic to S (semigroups 
with the same associated ideal) with nice system of generators. A system of 
generators is nice if it is formed by subsets of elements with disjoint support 
(see Remark [T^ . 

In section|3]decomposable varieties are introduced. These varieties represent 
geometrically the concept of decomposable semigroup. In this section we give a 
method to obtain a nice parametrization of a decomposable variety. 

In section [5l we show the relations between some generating sets of Is and 
, and we prove that Is is a complete intersection if and only if Is^ is a complete 
intersection for every i. Besides, we prove that a decomposable semigroup is a 
gluing if and only if at least one of the subsemigroups of its decomposition is a 
gluing. 

In section |6l we give a combinatorial characterization of the decomposable 
semigroups by using two different simplicial complexes, one of them defined in 
[2] and the other used in [TS] and introduced in [5]. 

In section [71 we illustrate our results with an example from Algebraic Statis- 
tics. In that example we reduce the number of indeterminates by using the 
decomposition of the semigroup. Furthermore we study the uniqueness of the 
Markov bases of the statistical models that appear in [M] and that have a de- 
composable semigroup associated. 

1 Preliminaries 

Let A = {ai, 02, . . . , a„} be the system of generators of S. Denote by I^a') 
the ideal of the semigroup {A') where A' is a subset of A. This ideal is in the 
polynomial ring in the subset of indeterminates of {Xi, . . . , X„} given by the 
elements of A' . Denote by k[^'] this subring of k[Xi, . . . , Xn] and its monomials 
hyXA'. 

Let supp (a) be the set {i\ ai ^ 0} where a G Z". In the same way, supp (X") 
is supp (a). For an element a G Z", a+ G Z" and G Z" are the unique 
vectors with disjoint supports such that a — — . These integer vectors 
have associated the binomial X°' — X°' G k[Xi, . . . , X„] whose support is 
defined as the union of the supports of and X" . 

In order to fix notation, for a decomposable semigroup S — {A) , the decom- 
position of S is the unique irreducible decomposition (up to permutations) 

5 (Ai) ® (A2) ® • • • ® (At) (2) 

where is the direct sum of subsemigroups and Ai is a set formed by some 
elements of A. If Si is the subsemigroup generated by {Ai) for i = 1, . . . , t, we 
obtain Is= hi^ ^ Ist- 
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For a lattice £ C Z", define its associated ideal as the binomial ideal gener- 
ated by - X"" \aeC\. 



The lattice of integer solutions of the linear system Ax = is denoted by 
keriS* (we identify the set A with the matrix (ai| • • • |a„)). The binomial ideal 
^kcrS C h[Xi, . . . , Xn] generated by the binomials 



is known as the lattice ideal of ker S and it is equal to the ideal Is (see [26l 
Lemma 9]). 

Consider the decomposition of ^ . The set ker S is the same lattice as 
ker((Ai U ^2 U • • • U At)) up to a permutation of some of its coordinates. To 
show this equality (up to permutations) we use the symbol = (for instance 
kerS* = ker((^i U A2 U • • • U At))). 

To study the uniqueness of Markov bases we need to study indispensable 
binomials (see jTj, [1], [TO], [S] and [TT]). The indispensable binomials are the 
binomials that belong (up to a scalar multiple) to every system of binomial 
generators of an ideal. There exists a unique Markov basis (up to a scalar 
multiple of its elements) of an ideal if and only if such ideal is generated by its 
indispensable binomials. In any case, to check if a Markov basis is unique it is 
necessary to perform calculations of high complexity in n indeterminates. 

Given the previous definitions, the ideal of a semigroup Is is a complete 
intersection if it is generated by n — rank (ker S) binomials (see [5]). The study 
of these kind of ideals is classical in commutative algebra and algebraic geometry. 

Other semigroups appearing in the literature are the semigroups that are 
the gluing of two. Gluing of semigroups was introduced by Rosales in [50]. The 
concept of gluing is an interesting tool to construct semigroups and ideals with 
some properties from semigroups and ideals with the same properties. Such 
properties are, for example, complete intersection or uniquely presented ideals. 

A semigroup S minimally generatecfl by Ci U C2 is the gluing of S' = (Ci) 
and S" = (C2), if there exists a set of generators of Is in the form 



where pi , p2 are sets of generators of Is' and Is" , and — X^ G Is with 
X" e k[Ci] and X^ G k[C2]. Equivalently, there exists d € S' f^ S" \ {0} 
such that G{S') n G[S") = G({d}), where G(S'), G{S") and G{{d}) are the 
associated commutative groups of S", S" and {d} (see Theorem 1.4 in [20] for 
details). Note that d is the S'-degree of X" . This element d G S* is called the 
gluing degree associated to the partition, in this case X" — X^ is called a gluing 



2 HNF-diagonalization of matrices 

In this Section, we use the Hermite normal form for matrices with integer coef- 
ficients. See details and algorithms in 5 . 

^We consider a minimal generator set of S because, in tlie otlier case, S is trivially the 
gluing of the semigroup generated by one of its non minimal generators and the semigroup 
generated by the others. 




piUp2U{X"-X''}, 



binomial. 
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Given a matrix B G TP^'^ ^ denote by F{Ii) C the set of rows of _B, and 
rF(B) = rowspanz(F(B)). 

Definition 1. We will say that an m x n matrix M = {rriij) with integer 
coefficients is in Hermite normal form ( abbreviated HNF) if there exists r < m 
and a strictly increasing map f from [l,m — r] to [l,n] satisfying the following 
properties: 

1. for I < i < m - r, ■mif(i) < I, mu ^ if I < i and mif(i) < rnkf(i) if 
k < i, 

2. the last r rows of M are equal to 0. 

By GLn (Z) we denote the group of matrices with integer coefficients which 
are invertible, i.e. whose determinant is equal to ±1. 

Theorem 2. Let L be an m x n matrix with coefficients in Z. Then there 
exists a unique m x n matrix H ~ (hij) in HNF of the form UL — H with 
U G GLn{Z). In this case, we write HNF(L)=H. 

For a given positive integer n and ^ < i, j < n,{i ^ j), we consider Cii^j, 
the matrix obtained from the identity matrix /„ interchanging the columns i 
and j. These matrices are known as column elementary matrices. We have 
det{Ci^^j) = — 1 and C^^j = C^J^j. Let M be an m x n matrix with integer 
coefficients. The matrix MCi^^j is the matrix resulting from interchanging the 
columns i and j in M . In the rest of this work, we denote this matrix as Q. 
We call this matrix, column permutation matrix. Similarly, we can define row 
permutation matrix. Observe that permutation matrices are nonsingular. 

Definition 3. We will say that an m x n matrix D with integer coefficients is 
an HNF-diagonal matrix if it satisfies the following conditions: 

1. the null rows are on the bottom side of the matrix, 

2. the null columns are on the right side of the matrix, 

3. every block Di has maximal rank, it is in HNF and for each disjoint par- 
tition of F{Di), Bi U i?2, we have: 

I y supp(/)]fi[ y supp(/)]/0. 

That is, an HNF-diagonal matrix has the following shape: 



( ^1 


e 


e ■ 


■ e 


\ 


e 




e ■ 


• e 




e 


e 


. 


• e 




V e 


e 


e • 




I 



(3) 



where for i — 1, . . . ,t — 1, Di is an HNF-matrix with non-zero columns, Q are 
the null matrices with the corresponding orders and Dt can be a null matrix or 
an HNF-matrix with zero rows and/or column^. 

^To compute it, we only consider the non-zero rows and columns of Dt- 
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For example, the matrix 



/ 1 





-1 


2 








\ 





4 


4 


























3 


1 


2 























V 

















/ 



is an HNF-diagonal matrix, where 

4 o) ^"'^ D,^{3 1 2). 

Definition 4. We will say that an m x n matrix L with integer coefficients 
is HNF-diagonalizable if there exists an unimodular matrix P and a column 
permutation matrix Q such that PLQ is an HNF-diagonal matrix. 

This is equivalent to saying, an m x ?i matrix L with integer coefficients is 
HNF-diagonalizable if there exists a column permutation matrix Q such that 
HNF(LQ) is an HNF-diagonal matrix. 



We have the following Lemma: 

Lemma 5. Let L be an m x n matrix with integer coefficients. Then L is HNF- 
diagonalizable if and only if there exists a column permutation matrix Q and an 
HNF-diagonal matrix D such that 

rF{LQ)^rF{D), 

and then rF{L) = rF(LQ) = rF(D) = tF{Di) x • • ■ x rF(A)- 

Proof. The proof follows from Test for Equality (see 5, P-74]) □ 

Thus a matrix L is HNF-diagonalizable if and only if the lattice rF(L) is, 
up to a permutation of its coordinates, the Cartesian product of some of its 
sublattices. 

Lemma 6. There exists an algorithm to determine if a matrix L is HNF- 
diagonalizable. If L is HNF-diagonalizable this algorithm finds two matrices Q 
and P such that PLQ is an HNF-diagonal matrix. 

Proof. By LemmaO L is HNF-diagonalizable if and only if there exists a column 
permutation matrix Q and an HNF-diagonal matrix D such that rF{LQ) = 
iF{D). Since Q is a column permutation matrix = Q and thus rF(L) = 
rF(DQ-i) = rF(DQ). Therefore HNF(L) = HNF(i:>Q). Trivially HNF(i:>g) is 
equal to DQ up to a permutation of its rows. 

In any case, the above condition can be detected by checking the supports 
of the rows of HNF(L). □ 

We can now describe an algorithm satisfying the previous Lemma. 

Algorithm 7. HNF-diagonalization 

In: An m x n matrix L with integer coefficients. 

Out: This algorithm detects if L is HNF-diagonalizable and finds two matrices 
Q and P satisfying Definition 
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1. We calculate H, the Hermite normal form of L, 



H 



( h,\ 

h2 



Let P2 be a unimodular matrix such that P2L = H. 

2. Initialize A — {2, . . . , m}, Aj is the j— th element of A and Bi — supp/i^, 
Fi = 9 for every i = l,...,m. Fi = {l} and i = 1. 

3. For j ^ 1 to card(A) (the cardinality of A) do: 

• If n Ba^ ^ 0, = 5, U Ba^, F,^F,U Aj and A = A \ A^. 

• J = J + 1- 

• If A = and i = 1 then L is not HNF-diagonalizable and terminate 
the algorithm. 

• If A = and i 7^ 1 then L is HNF-diagonalizable and go to Step S) 



Otherwise, i 



1 and go to Step [3] 



4. Let Pi be the row permutation matrix to place the rows of H according 
to the sets Fj, 

^-Pl.cardCFi) 
/lF2,l 



M = PiH 



!'2,card(F2) 



e U' 



) 



where Fj^t is the th element in Fj, and let Q be the column permutation 
matrix such that MQ is the HNF-diagonal matrix. Observe that the 
matrix Q is determined by the sets Bj modified in Step 3 and the zero- 
columns. 

5. The matrices Q and P = P1P2 are the matrices that we want to find. 

From now on, given a matrix L, we suppose that its associated HNF-diagonal 
matrix has the same form as ([3]). 
Example 8. Let L be the matrix 



2 / 



-1 


-4 


12 


4 


-3 


-2 


8 


-4 


-4 


-9 


-3 





-8 


-6 


7 


8 


12 


4 


1 


14 


8 


7 


16 


6 


2 


9 


14 


4 


-1 


-4 


-3 


-1 


-3 


-2 


-2 
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We have 







/ 1 











-1 


2 


\ 









4 








4 








HNF(L) = 










3 


1 








2 































^0 

















J 


By using Algorithm [71 we obtain that M = 


H, and the 




/ 1 
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1 
































1 
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1 















V 

















1 J 





Note that Q changes columns 3 and 4 for columns 5 and 6 in M. 
Hence, 



HNF 
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1 

















^ 




-1 


-4 


12 
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-3 


-2 


8 


\ 









1 



















-4 


-4 


-9 




3 





-8 


-6 




















1 










7 


8 


12 


4 


1 


14 


8 























1 







7 


16 


6 


2 


9 


14 


4 














1 














I 


-1 


-4 


-3 




1 


-3 


-2 


-2 


J 















1 





























I 




















1 ) 












/ 1 





-1 


2 





































4 


4 
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1 




2 





































































^ 



















J 















3 Decomposable semigroups 

In this section we will characterize the decomposable semigroups by means of the 
HNF-diagonalization of a matrix with rows are a base of the lattice ker S. Firstly, 
we present a technical lemma that will be a useful tool for this characterization. 

Lemma 9. Let S = (A) be a decomposable semigroup as in (0). The following 
are equivalent: 

1- Is = Is^ +Is2+---+Is,- 

2. ker S = ker 5*1 x ker ^2 x • • • x ker St ■ 

Proof. Assume that Is — Isi + + ■ ■ ■ + Ist ^-nd consider the lattices generated 
by the exponents of the binomials in a generator system of each ideal /g. G k[Ai] . 
By Chapter 12], for each i, the above lattice is equal to kerS'^. Therefore, 
ker S = ker Si x ker S'2 x • ■ • x ker St- 

Conversely, let ker S = ker x ker S'2 x ■ • • x ker St ■ Trivially, Is^ + Is^ + 
• ■ • + C Is- 
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Now we prove the other inclusion. Given — X^^ £ Is with supp (a) n 
supp (/3) = 0, then a — /3 G ker S. Since ker S = ker 5*1 x ker 5*2 x • • • x ker St, 
there exists some Ui £ ker Si with i — 1, . . . ,t, such that a — /3 = (ui, . . . , Ut). 
In this case, we deduce that a = {u'^ , . . . , uf) and /3 = (u^, . . . , u^). 

We denote -.^ uf , (3i :— , and then 

x^-x^^ xi\xx ---xz- • • • (4) 

where Xl] - X^] e Is, C Is tor i ^ 1, ... , t. 
Realizing operations in the binomial 

^Ai^A2 ■ ■ • ^At ~ ^Ai^A2 ■ ■ ■ ^At - 

- \^Ai ^ ^Ai)^A2 "'^At + ^aAa2 "'^At ~ ^Ai^A2 ' ' ' ^ At - 

- /' v"i v^*! v"2 V"* _L v^i /' v"2 V"* v^M Isl 

- l^Ai ~-^aJ^A2 "'-^At +^Ail^A2 ■■•"^At ~^A2 ■■■^AtJ -U 

= TO-x^;)Mi+x^; ((X- -x^^jx- ...X- ...x^: =a 

= {Xl\ -X^; )Mi+X^; ((X^^ - X%)M2 + (X^^ • ■ • - X% . . . X^j) = . . . 
= ^ A/; (x^^: - x^; ) e /^^ + + . . . + , with Afj e k[Xi , . . . , x„] . 

□ 

The aim of the following result is the characterization of decomposable semi- 
groups by using the HNF-diagonalization. 

Theorem 10. Let S be a semigroup and L be a matrix such that ker 5 = rF(L). 
Then S is a decomposable semigroup if and only if L is an HNF-diagonalizable 
matrix. 

Proof. Suppose that S" is a decomposable semigroup. By Lemma [9l 

ker S = ker 5*1 x ker 5*2 x • • • x ker St . 

Consider Li a base of each kerS'i, where Li is in HNF. 
Let L' be the matrix 



L' = 



( ^1 


e 


e • 


• e \ 


e 




e ■ 


• e 


e 


e 


La • 


• e 


e 


e 


e • 


• Lt 


V e 


e 


e • 


■ e ; 



Since L an L' span equivalent lattices, we have that there exists a column 
permutation matrix Q such that rF{LQ) = rF(L') and HNF(LQ) = HNF(L') 

^M, = X-l---XZ 
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(by Lemma O. It is easy to see that HNF(L') is an HNF-diagonal matrix, 
therefore L is HNF-diagonahzable. 

Conversely, suppose that L is HNF-diagonahzable. In this case, if D is the 
HNF-diagonal matrix associated to LQ, 

rF(LQ) = rF(i:>) = rF(i:>i) x rF(i:i2) x • • • x rF(A), 

where Di represents the blocks of the matrix D. Hence, 

Is — IrF(LQ) ~ A-F(_Di) + Ii-F(D2) + ' ' ' + IiF{Dt)- 

We conclude that S is decomposable and the subsemigroups Si of its ir- 
reducible decomposition are the semigroups generated by the column of A in- 
dicated by the elements of the corresponding blocks Di. In this case, Is^ = 

Ii-F{Di)- □ 

In view of the proof, we can say that given an HNF-diagonalization of L, 
we have the irreducible decomposition S = (Ai) © {A2) © • • • © (At), where for 
every j, Aj is the set formed by the elements of A corresponding to the blocks 
Dj of the associated HNF-diagonal matrix. 

Remark 11. Note that if the semigroup S is torsion free we can study the 
decomposition of the semigroup S applying the algorithm [7] to the matrix whose 
columns are formed by a system of generators of 5. In this case we obtain the 
decomposition operating directly in S. 

The examples of Section |4j Section [7] and Table [1] have been done applying 
directly the decomposition over their corresponding semigroups due to these 
semigroups are torsion free. 

Remark 12. Algorithm [7] gives us not only a method to obtain the minimal 
decomposition of a semigroup S (in case this decomposition exists), it can be 
also saw as a algorithm to obtain a nice system generators of a semigroup 
isomorphic to S. 

It we have a torsion free semigroup, we obtain a nice system of generators. If 
the semigroup is not torsion free, then we can obtain a nice system of generators 
from the lattices ker S by applying [211 Chapter 2] . 

4 Decomposable varieties 

It is well known (see [531 Chapter 4]) that every torsion free semigroup S = 
(ai, . . . ,a„) C has an afHne variety determined by the set of zeros of its 
ideal Is C k[Xi, . . . , X„]. We denote that variety by V{Is) parametrized by the 
equations Xi = with i = 1, . . . ,n and / G (k*)"*. 

If we consider the decomposable semigroup ([5]), we obtain that 

viis) = v{is,) n y(/sj n • • • n f(/sJ, 

where all the varieties V{Isi) C k""*^'"^*) are embedded in k". We say that a 
variety is a decomposable variety if it is obtained from a decomposable semi- 
group. 

Clearly the above parametrization depends on the generators of the semi- 
group S, thought the variety V{Is) is determined by the ideal Is (the variety 
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depends on the relations of the generators of S). Thus, systems of generators 
of isomorphic semigroups represent the same variety. 

From Remark [1^ wc obtain not only a decomposition when such decompo- 
sition exists, we can determine a nice reparametrization . 

Example 13. Consider the toric variety determined by the parametrization 



Xi = 


7-1 

n 


'2 




ll 


X2 = 


7-4 

n 


7-4 

'2 


11 


If 


2^3 = 


7l2 

n 


,~9 
'2 


712 

'3 


'4 


X4 = 


;4 
n 




'3 


'4 


X5 = 


7-3 

n 




^3 


'4 


xe = 


7-2 

n 




714 

'3 


714 

'4 


Xj = 


;8 
n 


j-6 


'3 


'4 



7-1 
'5 

7-4 
7-3 

'5 

7-1 

'5 

7-3 

'5 

7-2 

'5 

7-2 

'5 



and let C the semigroup associated to the parametrization generated by 
the columns of the matrix 



/ 


-1 


-4 


12 


4 


-3 


-2 


8 


\ 




-4 


-4 


-9 


-3 





-8 


-6 






7 


8 


12 


4 


1 


14 


8 






7 


16 


6 


2 


9 


14 


4 




V 


-1 


-4 


-3 


-1 


-3 


-2 


-2 


/ 



Now by applying Algorithm [7] and Remark[TTl we obtain that the semigroup 
S is the direct sum of 



/ -1 

-4 



-4 -3 -2 \ 



V 



16 

-4 



14 
14 



0' 



-2 / 



/ 12 

-9 
12 

6 

V -3 



-2 / 



and that it is isomorphic (it has the same ideal up to permutation of the vari- 
ables) to the direct sum of the semigroups Si and 6*2, 



fl 

4 





\ 



-1 
4 






2 \ 




/ 







/ \ 



3 1 2 



\ y 



Using this decomposition we can give nice reparametrization of the affine 



^Noto that the matrix A is the same that the matrix of Example [S] but in this example its 
columms are the generators of the semigroup S. 
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toric variety V, 

' xi = qi 
X3 = qi 

< X4 ^ 93 ■ 

X5 = qi^ qt 

x& = qf 
. 2:7 = ql 

Furthermore the ideal Is is the direct sum of the ideals Is-^ C C[xi, X2, x^, xq] 
and C C[x3,X4,xj], 

Is = (xl - Xe, X1X5 - X2) + {xl - xj, xl ~ X3) , 

and 

V = V {{xl - Xe, X1X5 - X2)) n V [{xl - XT, xl - xz)) . 

5 Some applications of the decomposition 

By Lemma [9l we can characterize some properties of ideal Is using the decom- 
position Is = Isi + -^52 + ■ • ■ + Ist induced by the decomposition of semigroup 
S. From now on, we will consider the notation fixed in 

Proposition 14. Let Is = Isi + Is2 + • ■ • + ISt ^ decomposition of Is- Then 
Is is generated by the disjoint union of the generators of each Is^ ■ 

Proof. Trivial, because the generator system of each /s. do not have common 
variables. □ 

Remark 15. The above result is also true for: Markov bases, universal Markov 
bases, Grobner bases, Graver bases and universal Grobner bases, and Betti{S) — 
ULi Betti{Si). 

With the previous remark, we obtain a corollary that allows us to study the 
uniqueness of Markov basis of Is using the Markov bases of the ideals Ig. . 

Corollary 16. Let Is — Isi + Is2 + ■ ■ ■ + Ist be a decomposition of Is- The 
following are equivalent: 

1- Is is generated by indispensable binomials, that is, Is has a unique Markov 
basis (up to scalar multiple of its elements). 

2- Isi is generated by indispensable binomials, for all i = 1, . . .t- 

Proof- Trivial by using Remark 1151 □ 

The following proposition allows us to determine if the ideal Is is a complete 
intersection by using the ideals /s; • 

Proposition 17. Let Is = Isi + IS2 + ■ ■ ■ + Ist be a decomposition of Is- Is 
is a complete intersection if and only if Is^ is a complete intersection for every 
i = l,-.-t- 



12 



Proof. By Theorem llOl wc have rankL^ — rankD.; and ranki = X]i=i ranki^ = 
Si=i rankDi. Additionally, n — ranki is a lower bound of the number of gener- 
ators oi Is, and card(Ai) — rankL^ is a lower bound of the number of generators 
of 75. , for i = 1 . . . , t. 

Suppose that Is is a complete intersection. In this case the Markov basis 
of Is has n — rankL minimal generators (i.e. lower bound is reached). Since 
n — rankL = X]*=i(<^^^'^(^i) ~ rankL^), the Markov basis of each 15. cannot 
have more than card(Ai) — rankL^ elements (see Proposition [T4t . Therefore Isi 
is a complete intersection for i = 1, . . . ,t. 

Conversely, if the Markov bases of the ideals /g. has card(^i) — rankLj 
elements, analogously the Markov basis of Is has n — rankL elements. Hence, 
Is is a complete intersection. 

□ 

The following result establishes that the gluing property can be characterized 
by the decomposition of the semigroup. 

Proposition 18. S is a gluing if and only if there exists i S {1, ■ ■ ■ ,t} such 
that Si is a gluing of semigroups. 

Proof. Suppose that S is the gluing of S' and S". In this case, there exists a 
gluing binomial X°' — e Is verifying that its S'-degree d G 5 is in 5" fl S" . 
Moreover, by 

where - X^; e Is, Cis ior i = 1, . . . , t. 

If we denote di & Si, the S'-degree of X"^', we have d = c?i H \-dt. Without 

loss of generality, we will suppose that di ^ 0. 

Clearly, di G {Si H S') n {Si H S") and di e G{Si n S") n G{Si n 5") c 
G{S') n G{S") = dZ. Since S n (-5) = {0} and d di + ■ • ■ + d*, we have 
di = d. Therefore, di = for every i > 1, and d = di ^ Si. That is, the gluing 
degree associated to the partition S = S' + S" belongs to a semigroup in the 
decomposition of S. 

As G{Si n S') n G(Si n S") = diZ and di e (Si n S') n (Si n S"), the 
semigroup Si is the gluing of Si fl S' and Si fl S". 

Conversely, if any Si is the gluing of S'^ and S'/ , then, it is easy to prove that 

S is the gluing of S,' and Sf U (u^,^, sA. □ 



6 Combinatorial results 

Following the idea of the introduction, it is clear that given a Markov basis B 
of an ideal of a semigroup Is such that it can be decomposed into two proper 
and disjoint subsets Bi and B2 such that B = BiU B2 and 

U supp(/) fl y supp(/) 

yfeBi J \feB2 

then the semigroup S is decomposable. A not necessarily irreducible decom- 
position is determined by the generators of the semigroup given by the sets 
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U/eBi supp (/) and U/ee, supp (/). The irreducible decomposition can be ob- 
tained applying this idea to the subsets Bi and 62- We have just seen, the 
decomposition of a semigroup is characterized by the decomposition into sub- 
sets with disjoint support of a Markov basis of its associated ideal. 

However, Markov bases can be characterized in a combinatorial way through 
the study of some simplicial complexes, such as the complex defined in [5] 
or the complex Vm used in |15) . where 

A™ = {f C {1, . . . , n}| m - ^ a, e 5}, 

V„i ~ |f C {monomials in k[Xi, . . . ,Xn] of ^-degree m}\ gcd(F) =/= l| 

and gcd(F) denotes the greatest common divisor of the monomials in F. 

A Markov basis of the semigroup ideal has an element of S'-degree to e S* 
if, and only if. A™ (respectively Vm) are not connected. Furthermore, the 
elements of a given degree, to G Betti{S), that appear in a Markov basis are 
fixed by the connected components of such complexes and the binomials of this 
degree can be obtained from Vm. For every non-connected complex Vm, the 
number of binomials obtained is equal to the number of connected components 
minus one (see [5] and [TS] for further details). 

Therefore, we have the following combinatorial characterization of decom- 
posable semigroups. 

Theorem 19. Let S be a semigroup. The following assertions are equivalent: 

1. S is decomposable. 

2. There exist Ci, C2 proper and disjoint subsets of {ai, . . . , a„} fulfilling: 

• {ai,...,a„} ^ Ci UC2. 

• For all m £ S with V„i not connected, all the vertices of V„i belong 
to k[Ci] or k[C2], 

• There exist toi,TO2 £ S with Vmi cind Vm^ '^oi connected such that 
the vertices o/Vm^ belong to k[Ci] and the vertices o/V„i, belong to 
k[C2]. 

3. There exist Ci, C2 proper and disjoint subsets of {1, . . . , n} satisfying: 

• {l,...,n} = CiUC2. 

• For all m £ S with A^ not connected, all vertices of A^ belong to 
Ci or C2, 

• There exist mi,TO2 £ S with A^^ and l^m2 '^oi connected such that 
the vertices of Ami belong to Ci and the vertices of A^^ belong to 

C2. 

Proof. Since the facets of A^ are the union of the supports of the monomials 
of the facets of Vm, ^ and 3 are equivalent. 

We now prove that 1 implies 2. Assume that S is decomposable. In such 
case, a Markov basis of Is is the union of the Markov bases of the ideals Jg. . 
With the notation fixed in consider Ci = Ai and C2 = Ui=2^j- ^ 
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be the S'-degree of any binomial of a Markov basis of any /g. (without lost of 
generality we consider z = 1). In this case Vm is non-connected and it has 
at least two different connected components, each of them with a monomial in 

If there exists Xg^X^.^ e in V™, then there exists Xg^X^.^ -^Ci ^ 

Is with S'-degree m. Therefore the element (a — 7, /3) belong to kerS"! x • • • x 
kerS't. Hence a — 7 S keriSi and /3 g kerS'2 x • • • x ker5t, that is not possible 
since S f^ {—S) = {0}. This concludes 1 implies 2. 

We now see that 2 implies 1. Consider the sets Ci and C2 given in the 
statement. In this case, since the Markov bases of Is are calculated from the 
connected components of the complexes Vm and such complexes are contained 
in k[Ci] or k[C2], these bases are formed exclusively by binomials in k[Ci] and 
k[C2]. By hypotheses, these bases have elements of both types. 

Clearly Is = I{Ci) + ^(02) ^"^^ ^ Grobner basis of Is is formed by the union 
of one Grobner basis of /(Ci) ^-nd another of I{C2)- 

It is straightforward to prove that S — (Ci) + (C2). It only remains to prove 
that this sum is a direct sum. If not, there would exist m' 6 (Ci) fl (C2) and 
a binomial X^_^ — X^^ £ Is with S'-degree m'. Since the normal form of this 
binomial with respect to a Grobner basis of Is is not equal to zero, this is not 
possible. □ 

Although from the above characterization one can obtain an algorithm to de- 
tect whether a semigroup is decomposable, this algorithm would not be efficient. 
We use combinatorics only to explain why the semigroup can be decomposed 
and not for computing the decomposition. 



7 Example "The binary graphical model of the 
bipyramid graph" 

We apply the results of this work to the statistical model The binary graphical 
model of the bipyramid graph (example BPg_bin in [14)). The semigroup of this 
model is determined by the columns of the matrix 



00000000000100000000000 
00000000000010000000000 



00000010000000100000000 
00000001000000010000000 



V 
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Note that the ideal I^a) belongs to a polynomial ring in 32 indeterminates. 
Applying Algorithm [7] to the matrix A (see Remark [TT|) , we get A is HNF- 
diagonalizable and its associated HNF-diagonal matrix is 



D = 



where Di = 



( 

e 
e 
e 
e 
e 
e 
e 

V e 

^10 
1 



e 

e 
e 
e 
e 
e 
e 
e 



e 


e 


e 


e 


e 


e 


e 


e 


e 


e 


e 


e 




e 


e 


e 


e 


e 


e 




e 


e 


e 


e 


e 


e 




e 


e 


e 


e 


e 


e 




e 


e 


e 


e 


e 


e 




e 


e 


e 


e 


e 


e 




e 


e 


e 


e 


e 


e 



-1\ 



and the permutation matrix Q swaps the 



\ 1 1 / 

columns of A grouping them in the following sets: 

Ai = {1, 9, 17, 25}, A2 = {2, 10, 18, 26}, A3 = {3, 11, 19, 27}, = {4, 12, 20, 28}, 

A = {5, 13, 21, 29}, Aq = {6, 14, 22, 30}, A7 = {7, 15, 23, 31}, = {8, 16, 24, 32}. 

In this example, we identify the index set Ai with the columns of A. 

Through an extremely fast calculation (which can be done manually), we 
obtain: 

• the Markov basis oi Ii^Ai) = I{Di) C k[a:i, xg, 0:17, X25] is {0:1X25 —xgxn}, 

• this basis is unique, 

• I{Ai) is a complete intersection, 

• (Di) is the gluing of the semigroups 5*^ = ((1, 0, 0), (— 1, 1, 1)) and 6*2 — 
((0, 1, 0), (0, 0, 1)), and then (Ai) is the gluing of the semigroups generated 
by the columns {1, 25} and {9, 17} of A. 

Therefore, by the results presented in this paper, the ideal I(^a) k[xi , . . . , X32] 
verifies: 

• has a unique Markov basis, 

{2:1X25 - X9X17, X2X26 - X10X18, X3X27 - a;iia;i9, X4X28 - X12X20, X5X29 - X13X21 
X6X30 - 2:14X22, X7X31 - X15X23, X8X32 - X16X24} , 

• it is a complete intersection. 

Furthermore, the semigroup S = (A) is the gluing of two semigroups and it 
decomposes as the direct sum of eight subsemigroups 

S = (Ai) © • • • e {As). 

Note that instead of performing highly complex computations in the ring 
of polynomials in 32 indeterminates, we performed the computations in rings 
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in four indeterminates. In the example and given the particular nature of the 
matrix D we need only to calculate a unique ideal. 

In Table m using jT^ Corollary 14], we complete the study on the uniqueness 
of the Markov bases for the examples that are decomposable in the database 
[141 . In addition, we collect the number of semigroups in the irreducible decom- 
position and the number of generators of each semigroup of the decomposition. 





N° of semigroups 


N° of generators 


Unique 


BM4r2-2_bin 


2 


8-8 


no 


G15g_bin 


2 


8-8 


YES 


G17g_bin 


4 


4-4-4-4 


YES 


5-4ml_bin 


2 


16-16 


yes 


BM5r3-2_bin 


2 


16 16 


no 


BM5r3-4_bin 


2 


16-16 


no 


BPg.bin 


8 


4-4-4-4-4-4-4-4 


YES 


G34g_bin 


2 


16-16 


NO 


G40g_bin 


2 


16-16 


NO 


G42g_bin 


2 


16-16 


YES 


G45g_bin 


2 


16-16 


YES 


G46g_bin 


4 


8-8-8-8 


no 


G47g_bin 


2 


16-16 


YES 


SPg.bin 


2 


16-16 


YES 


6-lI2_bin 


4 


16-16-16-16 


yes 


Glllg_bin 


2 


32-32 


NO 


G117g_bin 


2 


32-32 


NO 


G133g_bin 


2 


32-32 


NO 


G135g_bin 


2 


32-32 


NO 


G136g_bin 


2 


32-32 


NO 


G144g_bin 


2 


32-32 


NO 


G156g_bin 


2 


32-32 


NO 


G158g_bin 


2 


32-32 


NO 


G161g_bin 


4 


16-16-16-16 


NO 


G162g_bin 


2 


32-32 


NO 


G164g_bin 


2 


32-32 


NO 


G165g_bin 


2 


32-32 


YES 


G92g_bin 


2 


32-32 


NO 



Table 1: Decomposable semigroups of [Tl]. "no/yes" uniqueness solved in '14'. 
"NO/YES" uniqueness solved in this work and not in |14) . 



It is noteworthy that using the decomposition obtained in Table [TJ we low- 
ered the time for computing the Markov bases of the ideals by more that 50% . 
These times have been obtained using the program 4ti2 (see '27'). 

However, making the calculations with the library toric.lib of Singular 
(see [25| ) , the improvement in the computation of Markov bases was more than 
significant. For instance for G161_bin, using its decomposition the Markov basis 
was done in three seconds and without using the decomposition we stopped the 
computation after 24 hours with no result. 
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Remark 20. This paper can be applied to any finitely generated semigroup with 
zero element, cancellative (that is to say, if m + n = m + n' , with m,n,n' S 
S, then n = n') and such that 5" fl {—S) = {0}, since these semigroups are 
isomorphic to subsemigroups of an Abelian group (see Pl Section 1] for further 
details) . 
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